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Abstract. Although usual hierarchical clustering scenar- 
ios agree with the main observed properties of Lyman-a 
clouds and galaxies, they seem to fail to reproduce the 
cluster temperature - X-ray luminosity relation. A possi- 
ble explanation is that the IGM is reheated by supernovae 
or quasars before cluster formation. In this article, using 
a unified analytic model for quasars, galaxies, Lyman-a 
absorbers and the IGM, we obtain the redshift evolution 
of the temperature and the entropy of the gas and the 
corresponding cluster temperature - X-ray luminosity re- 
lation. We consider three scenarios: 1) no energy source in 
addition to photoionization heating, 2) heating from su- 
pernovae and 3) heating from quasars, in an open universe 
flm = 0.3, 17a = 0, with a CDM power-spectrum. We show 
that although quasars can easily reheat the IGM and raise 
its entropy up to the level required by current cluster ob- 
servations the energy provided by supernovae is unlikely to 
be sufficient. Indeed, the efficiency factor needed for the 
supernova scenario is of order unity {asN — 1-7) while 
for quasars we get aq ~ 0.008. Thus the IGM is more 
likely to have been reheated by quasars. Moreover, we find 
that if both scenarios are normalized to present observa- 
tions the reheating due to quasars occurs somewhat earlier 
{zs ~ 2) than for supernovae {zs ~ 0.4) because of the 
sharp drop at low z of the quasar luminosity function. We 
also show that the Compton parameter y induced by the 
IGM is well below the observed upper limit in all cases. 
Finally, we note that such a reheating process may partly 
account for the decline at low redshift of the comoving star 
formation rate and of the quasar luminosity function. In 
particular, we show that the contradictory requirements 
arising from clusters (which require a large reheating to 
affect the relation T — Lx) and galaxies (which require a 
small reheating so that galaxy and star formation are not 
too much inhibited) provide strong constraints on such 
models. Thus, the IGM should be reheated at low redshift 
z ^ 2 up to T ~ 5 10^ K. On the other hand, the reion- 
ization process of the universe is almost not modified by 
these entropy sources which means that our predictions 
for the former should be quite robust. 
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1. Introduction 

An important test of cosmological models is to check 
whether they can reproduce the wide variety of astro- 
physical objects we observe in the present universe. Sev- 
eral studies have already shown that the usual hierarchi- 
cal scenarios (such as the standard CDM model) provide 
predictions which agree reasonably well with observations 
for galaxies (Valageas & Schaeffer 1999a; Kauffmann et 
al.l993; Cole et al.l994), Lyman-a clouds (Valageas et 
al.l999a; Petitjean et al.l992; Miralda-Escude et al.l996; 
Riediger et al.l998), quasars and reionization constraints 
(Valageas & Silk 1999a,b; Gnedin & Ostriker 1997; 
Haiman & Loeb 1998). However, within this framework 
the simplest model leads for clusters to a temperature - X- 
ray luminosity relation Lx oc (Kaiser 1986) which dis- 
agrees with observations (Ponman et al.l996). It has been 
argued (Evrard & Henry 1991; Cavaliere et al.l997; Pon- 
man et al.l998) that this discrepancy could be explained 
by a "preheating" of the gas which would raise its entropy 
before clusters form. Indeed, this entropy "floor" would 
lead to a maximum density for the ICM which would break 
the previous self-similar scaling and provide a steeper re- 
lation T — Lx ■ On the other hand, it has also been sug- 
gested (Blanchard et al.l992; Blanchard & Prunet 1997) 
that the overcooling problem linked to galaxy formation 
(i.e. the fact that a small fraction of the baryonic content 
of the universe has been converted into stars while sim- 
ple estimates show that most of the baryons should have 
been able to cool by now) requires a reheating of the IGM 
in order to prevent the gas from cooling and falling into 
the dark matter potential wells. Thus, it is important to 
obtain a good handle on the "entropy history" of the uni- 
verse since it may play a key role in structure formation 
processes, for galaxies as well as for clusters. 
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In this article, we present an analytic model to derive 
the evolution of the entropy of the gas. To this order, we 
use a description developed in Valageas & Silk (1999a, b) 
to study the reheating and reionization history of the uni- 
verse by the radiation emitted by stars and quasars. This 
model also provides a consistent description of galaxies 
(Valageas & Schaeffer 1999a), Lyman-a clouds (Valageas 
et al. 1999a) and clusters (Valageas & Schaeffer 1999b). 
Thus, it ensures that we obtain a realistic scenario. Then, 
as detailed in Sect.H we add to this simple model an ad- 
ditional source of heating which corresponds to a direct 
energy input into the IGM from supernovae or quasars. 
We also describe the modifications induced by entropy 
considerations. Next, in Sect.^ we present the numerical 
results we obtain for an open universe for both super- 
nova and quasar heating. Finally, in Sect.^ we point out 
the feedback of such an entropy production onto structure 
formation, for galaxies, quasars and clusters. We also de- 
scribe the cluster temperature - X-ray luminosity relation 
we obtain from these scenarios. 

2. Model 

First, we briefly present the main characteristics of our 
model. Some of them are described in more details in 
Valageas & Silk (1999a) (hereafter VS). 

2.1. Multiplicity Junctions 

In order to evaluate the reheating and reionization of the 
universe by stars and quasars we need the mass func- 
tions of galaxies and QSOs. We also derive the multiplicity 
function of Lyman-a absorbers which provide most of the 
opacity at low z after reionization. This is necessary in 
order to model the density fluctuations in the IGM it- 
self. Indeed, at low z we consider density contrasts from 
(1-l-A) ~ 10"^ (voids) up to (H-A) 20 (massive Lyman- 
a forest clouds) within the IGM. Of course, the need to 
describe such a large class of objects, defined by various 
density thresholds (which will also be required by entropy 
considerations, as shown below) or other constraints (e.g. 
on their size), means that we cannot use the familiar Press- 
Schechter prescription (Press & Schechter 1974). Indeed, 
the latter is restricted to "just-virialized" halos defined 
by a density contrast Ac ^ 177. Thus, we assume that 
the non-linear density field is well described by the scal- 
ing model developed in Balian & Schaeffer (1989). This 
description is based on the assumption that the many- 
body correlation functions obey specific scaling laws which 
can also be seen as a consequence of the stable-clustering 
ansatz (Peebles 1980): 



^p(Ari,...,Arp) = A-^(^'-i)CpC 



-rp) 



(1) 



where 7 is the local slope of the two-point correlation func- 
tion. This model has been checked through the counts- 
in-cells statistics against various numerical simulations 



(e.g. Colombi et al.l997; Valageas et al. 1999b; Munshi et 
al.l999). Then, for a given class of objects defined by a re- 
lation A(M, z) (which implies a specific radius R{M,z)) 
we attach to each object the parameter x defined by: 



x{M, z) 



where 



1 + A(Af, z) 
l[R{M,z),z_ 



(2) 



r2 



y2 



6(i"i,i"2) with 



is the average of the two-body correlation function 
C2 (ri , r2) over a spherical cell of radius R and provides the 
measure of the density fluctuations in such a cell. Then, we 
write the multiplicity function of these objects (defined by 
the constraint A(M, z)) as (Valageas & Schaeffer 1997): 



viM,z)—^{-x^H{x) — 
Al M X 



(3) 



where p is the mean density of the universe at redshift z, 
while the mass fraction in halos of mass between M and 
M + dM is: 



fi{M,z)— ^x^H{x) — 
M X 



(4) 



The scaling function H(x) only depends on the initial 
spectrum of the density fluctuations and must be obtained 
from numerical simulations. In practice, we use the scaling 
function H{x) obtained by Bouchet et al.(1991) for a CDM 
universe. The mass functions obtained in this way have 
been checked against the results of numerical simulations 
for various constraints A(M) in the case of a critical uni- 
verse with an initial power-law power-spectrum (Valageas 
et al. 1999b). This description also takes into account the 
substructures which may exist within larger objects and 
it allows one to derive for instance the amplitude of the 
density fluctuations in the IGM (i.e. the "clumping fac- 
tor"). On the other hand, note that a simpler model where 
the density field is described as a collection of smooth ha- 
los with a universal density profile is inconsistent with 
the results of numerical simulations, as shown in Valageas 
(1999) (it implies a wrong behaviour of the many-body 
correlation functions). 

2.2. Galaxies 

Along the lines of VS, we use a simplified version of the 
model described in Valageas & Schaeffer (1999a). We de- 
fine galaxies by two constraints: I) a virialization condi- 
tion A > /S.c{z) (where Ac ^ 177 is given by the usual 
spherical model) and 2) a cooling condition tcooi < tn 
which states that the gas must have been able to cool 
within a few Hubble times tn at formation in order to 
fall into the dark matter potential well and form a galaxy 
(see also Rees & Ostriker 1977; Silk 1977). At late times. 
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the second condition 2) is the most restrictive for just- 
virialized halos (defined by A = Ac(z)) with a large 
virial temperature T . This means that these objects are 
groups or clusters which contain several subunits (i.e. 
galaxies) which satisfy the constraint 2) which approxi- 
mately translates in this case into a constant "cooling ra- 
dius" Rcooi ~ 100 kpc (see Valageas & Schaeffer 1999a 
for a detailed discussion). Thus, we define galaxies by 
the relation /S.gai{M,z) such that A^a; = ^c{z) except 
for high temperature halos at low z where this constraint 
would imply R > Rcooi- Then, these objects are defined 
by the condition R — Rcooi which gives their density con- 
trast Agai{M, z). This allows us to distinguish clusters or 
groups from galaxies. In practice, as seen in Valageas & 
Schaeffer (1999a), the condition 2) only plays a role at 
low redshift z < 1 for high temperature galaxies (T > 10^ 
K). Finally, we also require the virial temperature T of 
the galactic halos to be larger than a "cooling tempera- 
ture" Tcooi{z)- The latter corresponds to the smallest just- 
virialized objects which can cool efficiently at redshift z, 
defined by the constraint tcooi = tn- Moreover, it is con- 
strained to be larger than or equal to the temperature of 
the IGM. At low z we have Tcooi ~ 3 10'' K in our original 
model (VS) since for smaller temperatures cooling is very 
inefficient (due to recombination). From the lower bound 
Tcooi and the density threshold Agai{M, z) we obtain the 
galaxy mass function using (^). 

Then, we use a simple star formation model to derive 
the stellar content and the luminosity of these galaxies. 
This involves 4 components: (1) short lived stars which are 
recycled, (2) long lived stars which are not recycled, (3) 
a central gaseous component which is deplenished by star 
formation and ejection by supernovae winds, replenished 
by infall from (4) a diffuse gaseous component. The star 
formation rate dMs/dt is proportional to the mass of cen- 
tral gas with a time-scale set by the dynamical time. The 
mass of gas ejected by supernovae is proportional to the 
star formation rate and decreases for deep potential wells 
as 1/T, in a fashion similar to Kauffmann et al.(1993). 
Some predictions of this model (galaxy luminosity func- 
tion, TuUy-Fisher relation) have already been checked 
against observations (Valageas & Schaeffer 1999a). Thus, 
we obtain for the mean star- formation rate per Mpc'^: 



dps \ _ ^ p{z) 
dt I il„ tn 



X{x) e 



-X{x) 



x^H{x) — (5) 

X 



with: 



X{x) 



P_ 



1 



SN 



^ {\ + A)gai{x) 
(1 + Ae) 



(6) 



where p/Pd = 0.5 is a parameter of order unity which 
enters the definition of the dynamical time, while Tsn = 
10^ K describes the ejection of gas by supernovae and 
stellar winds (see also Kauffmann et al.l993): 

2 £SN EsN nrrip rjsN 



T: 



SN 



3 k msN 



10" K 



Here p is the mean molecular weight, csn ^ 0.1 is the 
fraction of the energy Esn delivered by supernovae trans- 
mitted to the gas {Esn = 10^^ erg) while rjsN/'msN — 
0.005 Mq^ is the number of supernovae per solar mass 
of stars formed (note that in VS we used Tsn — 2 10^ 
K). The factor A e^^ in (^) comes from the dependance 
of the efficiency of star formation on the properties of 
the host galaxy. Thus, small galaxies with a shallow po- 
tential well (T ^ Tsn) are strongly infiuenced by su- 
pernovae and stellar winds which eject part of the gas 
so that a small fraction of the baryonic matter is con- 
verted into stars (A <C 1). On the other hand, at low z 
large halos (T ^ Tsn) which formed at a high redshift 
{Agai > Ac) have already converted most of their gas into 
stars (A ^ 1). Indeed, their high density (due to their 
large redshift of formation) translates into a small dynam- 
ical time, hence to very efficient star formation within our 
model, which leads to the factor ^^(1 -I- A)gai in (||). As 
in the model used in Valageas & Schaeffer (1999a) and 
VS we assume that the gas ejected from the inner parts 
of small galaxies (T ^ Tsn) remains bound to (or close 
to) the galactic halo so that it can cool and fall back into 
the galaxy. This leads to a self-regulated star formation 
process (within each individual galaxy) which takes care 
of the "overcooling problem". Note that an alternative, 
as suggested in Blanchard et al.(1992) (also Prunet & 
Blanchard 1999), would be that this gas gets mixed with 
the IGM and leads to a progressive heating of the IGM 
which prevents most of the gas to cool and form galax- 
ies. This will also correspond to our supernova heating 



scenario (SN), see Sect. 2.5 



2.3. Quasars 

In a fashion similar to Efstathiou & Rees (1988) and 
Nusser & Silk (1993) we derive the quasar luminosity 
function from the multiplicity function of galactic halos. 
Thus, we assume that the quasar mass Mq is proportional 
to the mass of gas Mgc available in the inner parts of 
the galaxy: Mq = F Mgc. In our case this also implies 
that Mq ~ F Ms where Ms is the stellar mass. We use 
F = 0.008 which is consistent with observations (Magor- 
rian et al.l998 find that Mq ^ 0.006 Ms). We also assume 
that quasars shine at the Eddington limit so that their life- 
time is given by ig — 4.4 eg 10® yr where eg = 0.1 is the 
quasar radiative efficiency and that a fraction Aq =0.1 of 
galactic halos actually host a quasar. Thus, the luminosity 
of a quasar of mass Mq is: 



Ln = 



£Q Mq 
to 



(8) 



and the quasar multiplicity function i]Q{MQ)dMQ/MQ is 
obtained from the galaxy mass function r]g{M)dM/M by: 



(7) VqIMq) 



dMQ 
Mq 



Min 



tM 



Aq %iM) 



dM 
If 



(9) 
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Here Im is the evolution time-scale of galactic halos of 
mass M defined by: 



t 



1 



M 



d 



rig{M) dt 



(10) 



Since the quasar life-time tq ^ 10^ yr is quite short, we 
have r]Q{MQ)dMQ/MQ = Aq tg dr]g/dt dM/M. This also 
means that r]Q{MQ)dMQ/MQ - Aq tQ/tn Vg dM/M. 
The factor tq/tH shows that the quasar luminosity func- 
tion is biased towards large redshifts as compared with 
the galaxy luminosity function. In particular, it peaks at 
z ~ 2 and shows a significant drop at smaller redshift 
while the galaxy luminosity function keeps increasing un- 
til z ^ and the star formation rate only decreases after 
z ^ 1. As we shall see in Sect.^ this implies that quasar 
heating of the IGM occurs earlier than supernova heat- 
ing. Note that we only have two parameters: (eq F/tq) 
and (Aq tq). Hence a larger fraction of quasars Aq with 
a smaller life-time tq would give the same results. More- 
over, the assumption that quasars shine at the Eddington 
limit gives the ratio (eq/tq) while the parameter F is 
constrained by the observed ratio (quasar mass) / (stellar 
mass). The normalization factor (Aq tq) is constrained 
by the observed quasar luminosity function. Our results 
agree reasonably well with available B-band observations 
for 0.16 < z < 4.5 (VS and SectJT^). 



2.4- Lyman-a clouds 

We also include in our model a description of Lyman-a 
clouds. These correspond to density fluctuations in the 
IGM as well as to virialized halos which may or may not 
have cooled. More precisely, we consider three different 
classes of objects (Valageas et al. 1999a). 

Low-density mass condensations with a small virial 
temperature see their baryonic density fluctuations erased 
over a scale Rd as the gas is heated by the UV background 
radiation (or other processes) to a temperature T forest ■ 
More precisely, we define the scale Rd by: 




ikTi 



IGM 



fitrip 



(11) 



where Cs is the sound speed, Tjqm the IGM temperature, 
tn the age of the universe, mp the proton mass and 7 ^ 
5/3. These mass condensations form a first population of 
objects, defined by the scale i?^, which can be identified 
with the Lyman-a forest at low z. We set the characteristic 
temperature T forest by: 



Z < Zri : T forest = Max (3 10^K,T/GAf) 
Z > Zri : T forest — TjGM 



(12) 



where Zri is the reionization redshift and Tjqm is the tem- 
perature of the IGM. At low z the term 3 10^ K models 
photoionization heating for the clouds (the IGM is also 



heated by the UV flux but in addition it undergoes adia- 
batic cooling because of the expansion of the universe, so 
that at low redshift z < 1 we can have Tjgm < 3 10* K). 
As explained in Valageas et al. (1999a), these absorbers 
are not necessarily spherical clouds of radius Rd- Some 
may be long fllaments with a length L ^ Rd and a thick- 
ness Rd- Moreover, they can also be interpreted as den- 
sity fluctuations within the IGM rather than distinct en- 
tities. Next, potential wells with a larger virial tempera- 
ture T > Tforest do not see their baryonic density pro- 
file smoothed out. Thus, they define a second class of ab- 
sorbers which for T > Tcooi correspond to the galactic 
halos (objects with Tforest < T < Tcooh if such a range 
exists, simply are virialized objects which have not cooled 
hence have not formed stars). Note that one such object 
can produce a broad range of observed column densities 
depending on the impact parameter of the line of sight. 
This population corresponds to Lyman- limit systems. Fi- 
nally, the deep cores of these halos are neutral because of 
self-shielding and they form our third class of absorbers, 
corresponding to damped systems. 

2.5. Evolution of the IGM 
2.5.1. Temperature evolution 

As in VS the gas in the IGM is heated by the back- 
ground radiation while it cools because of the expansion 
of the universe and several atomic processes (collisional 
excitation, ionization, recombination, molecular hydrogen 
cooling, bremsstrahlung and Compton cooling or heating) . 
Meanwhile, hydrogen and helium are reionized by the UV 
flux. In our calculation, we take into account the opacity 
due to the gas present in the underdense regions which 
fill most of the volume as well as the absorption due to 
discrete clouds (the "Lyman-a clouds" described above). 
We also follow the evolution of the HI, Hell and Helll fill- 
ing factors describing the ionized bubbles around galaxies 
and quasars, as well as the clumping of the gas (which also 
enters explicitly into the model for Lyman-a clouds). The 
evolution of the background radiation field J^, is obtained 
from the radiation emitted by stars and quasars, which we 
described in the previous sections, see VS for details. We 
write the evolution of the temperature of the IGM as: 



dTiGM ^ rr TiGM 
IT = - ±IGM - — 

at a tcooi 



theat.J 



IGM 
tE 



(13) 



where a{t) is the scale factor (which enters the term de- 
scribing adiabatic cooling due to the expansion) . The heat- 
ing time-scale theat,j which corresponds to photoioniza- 
tion heating is given by: 



47r 



heat, J 



3/2nf)kTj 



IGM 



E 



dv 



(14) 



where j — (HI,HeI,HeII), Vj is the ionization threshold 
of the corresponding species, n,- its number density in 
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the IGM and Ub the baryon number density. The cool- 
ing time-scale tcooi describes coUisional excitation, coUi- 
sional ionization, recombination, molecular hydrogen cool- 
ing, bremsstrahlung and Compton cooling or heating. We 
compute the redshift evolution of the ionization state of 
hydrogen and helium and we use the cooling rates from 
Anninos et al.(1997). In particular, as shown in the upper 
panel of Fig. 4 in VS at low z after reionization the main 
cooling processes in the IGM are adiabatic and Compton 
cooling. Indeed, when the medium is reionized coUisional 
excitation cooling is strongly suppressed (see discussion in 
VS and Efstathiou 1992). 

Finally, we added to the evolution equation we used 
in VS a new term T/gm/^e- This corresponds to an addi- 
tional source of energy, which we assume here to be uni- 
form. In particular, this term models in our framework the 
energy output provided by supernovae or quasars, which 
has been advocated in the litterature in order to raise the 
entropy level of the IGM (e.g. Ponman et al.l998; Tozzi 
& Norman 1999). As explained in Sect. 2.2, in our orig- 



inal model the influence of supernovae was restricted to 
their parent galaxy (note that this is consistent with nu- 
merical simulations by Mac Low & Ferrara 1999 which 
suggest that gas ejection is negligible for galactic halos 
with M > IO^Mq). In contrast, one model we investigate 
in this article corresponds to a "maximally efficient" sce- 
nario where the energy produced by supernovae reheats 
the IGM as a whole. In the actual universe, the effect of 
supernovae is likely to lie somewhere in-between these two 
cases, but these two models allow us to get an estimate 
of the allowed range for the reheating process (see also 
Tegmark et al.l993 for a study of reheating and reioniza- 
tion of the IGM by supernovae-driven winds) . 

We note that using a uniform source of energy (i.e. 
we do not let the energy source term vary in space as a 
function of the distance to the nearest galaxy or quasar, al- 
though we model this effect for photoionization heating) is 
probably a better approximation than it may seem at first 
sight. Indeed, at late times 2^2 when this process domi- 
nates most of the matter is embedded within positive den- 
sity fluctuations (filaments, virialized halos, see VS) which 
show a strong clustering pattern as seen in numerical sim- 
ulations (e.g. Bond et al.l996). Note that this is included 
in our model of the density field, described in Sect. 2.1. For 



instance, we obtained in Valageas et al. (1999a) the am- 
plitude of the two-point correlation function of Lyman-a 
clouds and we describe in Valageas et al. (1999c) the bias of 
the various objects we observe in the universe (Lyman-a 
clouds, galaxies, quasars, clusters). Thus, most of the vol- 
ume consists of low-density regions while most of the mat- 
ter is embedded within small or thin structures (filaments, 
halos) which are located close to galaxies since most clus- 
ters and galaxies form on density peaks within these mass 
condensations, though there may also be some isolated 
galaxies amid low-density regions (note that this "bias" 
translates into the correlations of these objects). As a con- 



sequence, the energy provided by supernovae or quasars 
does not need to travel very far in order to heat most of 
the matter. Indeed, for this it is sufficient to "spread" the 
energy over filaments while leaving cool voids in between. 
In this case, the temperature Tjgm would rather corre- 
spond to a "mass-averaged" temperature, describing the 
network of halos and filaments which contain most of the 
matter while voids would be cooler. Of course, one may 
also expect voids to be easily heated to the temperature 
of the filaments since due to their low density and small 
mass they only require a small amount of energy in or- 
der to reach the temperature of the neighbouring regions. 
Thus, the assumption of a uniform energy source appears 
to be a reasonable first order approximation. However, it 
is clear that a careful! study of this problem would be in- 
teresting, but this would probably require very detailed 
numerical simulations which are beyond the scope of this 
study. 

In this article, we consider the additional energy de- 
scribed by the term TjcM/tE in (p^ ) to be provided by 
supernovae or quasars. Thus, we can write: 



Tig 



M 



IGM 



IGM 



tE 



tE,SN tE,C 



(15) 



which explicitly shows these two possible sources of en- 
ergy. Using our model for galaxies which we described in 
Sect. 2.2, we can write the source term Tigm /tE,SN due 
to supernovae as: 



Tigm oisn T; 



SN 



tE,SN 0.1 Pbiz) 



dt 



(16) 



where Pi,{z) = ilt/^m 7>{z) is the mean baryonic density of 
the universe (the fraction of matter within stars is always 
negligible) and asN is the efficiency factor, similar to csn 
in (^ , which measures the fraction of the energy produced 
by supernovae which is available to heat the gas. Thus, we 
have asN < 1- Next, from the model of quasars presented 



in Sect. 2.3 we have for the quasar contribution: 



T 



IGM 



1 



tE.Q 0.1 3/2nbk 



dLc 



(17) 



where ag < 1 is the efficiency factor similar to eg in (H). 
However, if there are some additional energy sources (e.g. 
the decay of some exotic particles) we could have an effec- 
tive a larger than unity. Of course, in this case the time- 
dependence of this hypothetic energy source is unlikely 
to be proportional to the star or quasar formation rate 
and one should explicitly detail the origin of this process 
to get its time-evolution. In this article, we shall restrict 
ourselves to the formulation (^5|) which models the possi- 
ble effect of supernovae or quasars on the IGM, but one 
cannot disregard the fact that our source term TiQM/tE 
may in fact correspond to some new process. From the 
expressions (^6|) and ( |l7| ) we can directly obtain a simple 
estimate of the magnitude of these effects. Indeed, from 
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(|l^ we see that supernovae heat the IGM to a tempera- 
ture TiGM SN oi the order: 



TlGM,SN 



0.1 



Fstar TsN ^ C^SN K 



(18) 



where we used (0) and the fact that at late times z < 1 the 
fraction of baryonic matter which has been converted into 
stars is Fstar 0.1. Thus, supernovae can reheat the IGM 
up to 10^ K at most. On the other hand, from ([l7|), (||) 
and (^ we see that quasars heat the IGM up to T/ga/.q 
of the order: 



^qF 



Q 



Fstar Tq ^ Uq 10« K(19) 



where we used the parameters introduced in Sect, 2.3 and 
we defined: 



T, 



Q 



3fc 



4.3 10^2 K 



(20) 



The factor Fstar comes from the fact that Mq = F Mgc ^ 
F Ms, which agrees with observations. Thus, quasars can 
potentially heat the IGM to a very high temperature, much 
larger than the temperature induced by supernova heat- 
ing, because quasars are very efficient engines to convert 
the rest mass energy of matter into radiation or energy 
while a small fraction of the matter converted into stars 
leads to supernovae (~ 10~^) which themselves have a 
small efficiency factor (~ 10~^) so that Tsn ^ IO^^Tq. 
Note that the estimates ( p^ ) and (|l^) are very robust, 
independently of the details of the model, since they are 
directly constrained by the observed galaxy and quasar 
luminosity functions. On the other hand, the new param- 
eters as AT and ag are only constrained to be smaller than 
unity. One would need a detailed study of many physi- 
cal processes which are still poorly known to set a precise 
value for these efficiency factors. In this article, we shall 
treat them as free parameters, which we take to be con- 
stant in time. Thus, our goal is to evaluate the possible 
effects of these energy sources, which in turn will give us 
some constraints on their magnitude. 



2.5.2. Entropy evolution 

From the model of the IGM described in the previous 
sections, we can also obtain the evolution of the entropy 
of the gas. The entropy of a Maxwell-Boltzmann gas is 
given by the Sackur- Tetrode equation: 



S 



M.B 



(N, V, T) = Nk 



In 



(21) 



where A'^ is the number of particles, within the volume V , 
and: 



At = 



'27r?i^ 



(22) 



Thus we define the specific entropy S as: 



5 = log 



kT 



-2/3 



1 keV cm 



(23) 



where is the baryonic number density (and we note 
log the decimal logarithm). As explained in details in VS, 
we consider that at late times most of the volume of the 
IGM consists of large underdense regions with a density 
contrast (1 -f A)„ ("u" for underdense) given by: 



(1 -t- A)„ = Min 1 , ^{Rd) 



(24) 



This simply states that at high z (when £,{Rd) ^ 1) we 
have p„ —'p (i.e. the universe is almost exactly a uniform 
medium on scale Rd) while at low z we have <'p since 
most of the matter is now within overdense objects (clus- 
ters, filaments, etc.) while most of the volume is formed 
by underdense regions. The scale Rd was defined in (|ll]) 
while the exponent w is given by the power-law behaviour 
of the scaling function II{x) at small x. In addition to 
these "voids" and the virialized halos which are identi- 
fied to galaxies or clusters, there are also density fluctu- 
ations (clouds, filaments) which form the Lyman-a forest 
or small virialized halos which have not cooled (T < Tcooi) 
and can have a larger temperature than the underdense 
regions (the temperature of the gas within these small ha- 
los is of the order of the virial temperature of the potential 
well). It is important to take into account these density 
fiuctuations since at late times {z < 1) the "overdensity" 
(1 -I- A)„ is as low as ~ 10~^ while the density contrast of 
forest clouds reaches ^ 20. This wide variation of the local 
density means that the average entropy of the IGM can 
be significantly different from the entropy which would be 
computed from (23) within "voids". Thus, we first define 
the entropy Su characteristic of the underdense regions 
which fill most of the volume by: 

-2/3 \ 



Su = log 



kT, 



IGM riu 



1 keV cm 



where is the baryonic density obtained from (C 
we define a "mean IGM entropy" {S)igm by: 



(25) 



Then, 



(S) 



IGM 



(1 + A)/GM 



(1 + A)„ Su 



(26) 



S{x) x'H{x)- 
JO ^ \ 

where the mean density contrast (1 + A) jgm is given by: 
(1 + A),GA/ = (1 + A)„+ / x'H{x)— (27) 



It corresponds to the mean density (total mass over the 
volume) of the matter which is not embedded within viri- 
alized halos which have cooled, from very underdense re- 
gions up to forest clouds seen as density fluctuations in the 



p. Valageas & J. Silk: The entropy history of the universe 



7 



IGM. Since the entropy is an additive quantity, the rele- 
vant quantity is indeed the mean {S)igm which describes 
the average entropy of the IGM, see (^). In particular, at 
late times the quantity Su which corresponds to the small 
fraction of matter located within voids is much larger. We 
also define the mass-averaged overdensity {1 + A)Ly char- 
acteristic of the matter outside galaxies and clusters by: 



(1 + A)l, 



(1 + A)^ 



(1 + A) IGM 



(28) 



The quantity {1 + A)Ly corresponds to the average over- 
density of IGM particles, weighted by the number of par- 
ticles and not by the volume they occupy. Finally we in- 
troduce the mean temperature T^j,: 



Ly 



(1 + A)„ TiGM 




This takes into account the fact that some of the gas out- 
side galaxies and clusters is located in Lyman-a forest 
clouds with T > 3 10'* K (due to photoionization heat- 
ing) which may be larger than Tjgm (which also involves 
adiabatic cooling) , and possibly within some virialized ha- 
los with TiGM < T < Tcooi which have not cooled (due 
to their low virial temperature which implies inefficient 
cooling). 

2.5.3. Effect of the IGM entropy on galaxy formation 

As gravitational clustering builds increasingly large struc- 
tures, the baryonic matter content of the universe gradu- 
ally becomes embedded into virialize d h alos where it cools 
and forms stars, as described in Sect. 2. 2 where we detailed 
our model for galaxy formation. In particular, the "cool- 
ing temperature" Tcooi{z) which characterizes the smallest 
virialized halos which can cool was given by the condition 
tcooi = tn where the cooling time which depends on the 
density and the temperature of the gas satisfies: 



icool ^ 



T 



nb^{T) 



(30) 



where A(r) is the cooling function. Thus cooling is more 
efficient for larger baryonic densities rn, (because collisions 
are more frequent). In the original model the cooling time 
attached to a given halo was computed using the virial 
temperature Tyir for T and a density contrast Ac{z) ^ 177 
to obtain the gas density. However, if the IGM is preheated 
and gets a large entropy at earlier times, the gas may not 



follow the dark matter to form a mass condensation with 
mean baryonic density = {1 + Ac)pf,. Indeed, during 
the adiabatic collapse of the gas (before it cools) its tem- 
perature increases as Tad oc n^^^ ■ Hence the compression 
will stop if Tad reaches Tyir before the density contrast 
reaches Ac. Indeed, gas with T > Tyir does not fall into 
the potential well. Thus, we obtain an upper bound nb,ad 
for the baryonic density reached within a virialized halo 
of temperature T„ir, defined by: 



IGM 



= log 



kT 



-2/3^ 



vir "'b,ad 

1 keV cm^ 



(31) 



We can also write the density contrast (1 + A)ad given by 
(ill) as: 

(l-|-A),, = (l-fA)„ (I^Y\oiis.~iS).oM) (32) 

V ^ IGM J 



Thus, in order to compute the cooling time tcooi from (30) 
we use the overdensity (1 + A)b given by: 



(29) Ab = Min(Ac,A„rf) 



(33) 



Using this prescription we compute the "cooling temper- 
ature" Tcooi{z). It is obvious from ( ^ ) and ( |33| ) that the 
effect of the entropy of the intergalactic gas is to make 
cooling less efficient, which leads to a possible increase of 
the characteristic temperature Tcooi- In particular, at late 
times (z < 3) if the entropy production is sufficiently large 
it may happen that the cooling temperature Tcooi does not 
exist any more. Indeed, as we explained above the large 
entropy of the gas diminishes the gas density which enters 
(^0|). Moreover, the ratio T/A(T) displays a minimum at 
a finite value T ^ 10^ K since at large temperatures the 
cooling function behaves as A(T) oc y/T (bremsstrahlung 
is the main cooling process) while below 10''' K it nearly 
goes to (note that in our original model the high tem- 
perature halos which cannot cool are identified to clus- 
ters while the low temperature objects are Lyman-a ab- 
sorbers). As a consequence, if the entropy level is such 
that the gas density within the halos with a virial tem- 
perature ~ 10^ K is too small to allow efficient cooling, 
no halo can cool. Of course, this does not mean that there 
are no galaxies ! It simply means that all just- virialized 
halos (i.e. overdensities defined by the density threshold 
Ac) are identified to "clusters" (or "groups") in the sense 
that they consist of one or several smaller higher-density 
subunits (galaxies), which could cool at earlier times when 
they formed, embedded within a larger structure contain- 
ing some hot gas. 

Let us note zs the largest redshift where no halo can 
cool (i.e. Tcooi{zs) as defined above does not exist). We 
shall have < 3 since the entropy production is linked 
to galaxy or quasar formation and cooling is less efficient 
at low redshift where the baryonic density is lower. More- 
over, since the mean entropy {S)jgm increases with time. 
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at smaller redshifts 2; < we have the same situation. 
Then, along the lines developed in Valageas & Schaeffer 
(1999a) to distinguish galaxies from clusters at z ~ 0, we 
assume that gravitational clustering is stable and that, 
to a first order approximation, galaxies which form at zs 
do not evolve significantly at later times when they get 
embedded within larger structures which cannot cool as 
a whole. Thus, after zs these galaxies may get closer to 
form a group but we neglect their possible mergings. Note 
that the gas which cooled before zs and fell into these po- 
tential wells to build these galaxies formed "small" dense 
entities (the baryonic distribution extends to smaller radii 
than the underlying dark matter halo) which are likely to 
keep their identity for a longer time than their surrounding 
dark matter halos which may join to make a larger object. 
Furthermore, we note that the presence of substructures 
within dark matter halos themselves and the dependance 
of the characteristic density of a halo on its mass (it is pro- 
portional to the average density of the universe at the time 
this mass-scale turned non-linear, e.g. Navarro et al.l996) 
suggest that this picture may also be a good approxima- 
tion for the dark matter density fluctuations themselves 
(see discussion in Valageas 1999). Thus, we assume that 
after zs these galaxies keep their mass, radius and density 
unchanged. As a consequence, their density contrast at a 
later time is not lS.c{z) but: 

z<zs : (1 + A),,,(^) = (l + Ac)(zs) (1^)' (34) 

Thus, at small redshifts 2: < we no longer define galax- 
ies by the virialization condition A = Ac(z). Instead, we 



use the constraint A = Ago/(z) as defined in 



As ex- 



plained in Sect.^Tj this can be done in a straightforward 
fashion within our description of the non-linear density 
field: we simply use this density contrast Agai(z) in (||) 
to obtain the multiplicity function. Of course, as can be 
checked in (||) and (Q) the parameter x attached to such 
galactic halos does not evolve with time which also im- 
plies that the fraction of matter embedded within these 
objects is constant with time. Thus our prescription is self- 
consistent. This relies on the fact that in the non-linear 
regime relevant for galaxies the two-point correlation func- 
tion grows as £,{R,t) cx a{tY at fixed physical length i?, 
where a{t) is the scale-factor, as predicted by the stable- 
clustering ansatz (Peebles 1980). This behaviour is indeed 
consistent with numerical simulations (e.g. Valageas et 
al. 1999b). Then, at low redshifts we define Tcooi as being 
equal to the value it had at zs (more exactly for z ), 
with the constraint that it is larger than Tjqm'- 



z < zs ■■ Tcoo/(z) = Max [rcoo;(zs),T/GAf] 



(35) 



This is indeed the virial temperature of the smallest galax- 
ies, which cooled at zg, with T > Tjgm- This latter con- 
straint is due to the fact that the gas within small halos 
with T < TjQj^i will be heated up to Tjqm and it will 



escape from the potential well (but of course there will 
remain a small galaxy made of old stars). However, this 
condition does not play any role in practice since Tjgm 
does not increase much after zs- 

3. Numerical results 

We can now use the model we described in the previous 
sections to obtain the entropy history of the universe, as 
well as a consistent description of its reheating and reion- 
ization, together with the formation of quasars, galaxies 
and Lyman-a clouds. We shall consider the case of an open 
universe flm = 0.3, Ha = 0, with a CDM power-spectrum 
(Davis et al.l985), normalized to erg = 0.77. We choose 
a baryonic density parameter Qf, — 0.03 and Hq ~ 60 
km/s/Mpc. 

We shall consider three cases for the "entropy sce- 
nario". First, for reference we present some results we 
obtain for asN — ctQ = 0, as in VS. Then, we study 
both cases where only one of these two efficiency factors 
is non-zero. This allows us to see clearly the influence of 
each of these processes, as well as the magnitude of the 
relevant parameter a needed to get an appreciable effect. 
More precisely, in both cases we choose the value of a such 
that the mean IGM entropy {S)igm defined in (26) which 
we obtain at z = satisfies {S)igm — 2. Indeed, Ponman 
et al.(1998) find that the "entropy" s = T/nl^^ of small 
cool clusters seems to depart from the expected scaling 
law and to converge towards a floor value s ~ 100h~^^^ 
keV cm^. This provides an upper bound for the IGM en- 
tropy {S)igm[z — 0) since these objects have just formed 
and we can expect the entropy of the hot gas to increase 
(e.g. through shocks) rather than decrease during gravi- 
tational collapse (before it cools). Moreover, if we assume 
that supernovae or quasars are indeed the source of this 
entropy floor (since gravitational collapse effects shoud 
not break the expected scaling law) this gives the value 
of the corresponding parameter a. In the actual universe 
it might happen that both sources of heating, supernovae 
and quasars, have the same magnitude. Then, the param- 
eters asN and aq would be close to those we obtained for 
the individual cases. However, such a coincidence would 
be somewhat surprising. 

3.1. Reheating of the universe 

We show in Fig.|l| the reheating history we obtain for the 
three cases. The temperature Tm is a mass-averaged tem- 
perature which takes into account all the matter, from 
"voids" up to filaments, galaxies and clusters. The de- 
crease with time of the IGM temperature at large redshift 
z > 24 (log(l + z) > 1.4) is due to the adiabatic expan- 
sion of the universe. Then, at z < 24 the IGM is slowly 
reheated by stars and quasars until it reaches at z '--^ 9 a 
maximum temperature Tjgm ~ 3 10'* K where collisional 
excitation cooling prevents any further increase. Eventu- 
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Fig. 1. The redshift evolution of the characteristic tem- 
peratures of the universe. We display the IGM tempera- 
ture TjGM (solid curve) and the virial temperature Tcooi 
(dashed line) of the smallest galaxies. The upper panel 
corresponds to asAr = ag = (only photoionization heat- 
ing) , the middle panel labelled "SN" to supernova heating 
{ctSN 7^ , ag = 0) and the lower panel labelled "QSO" 
to quasar heating {asN = , aq ^ 0). The vertical line 
in both lower figures shows the redshift zg. In the upper 
panel we also show the mass-averaged temperature 
(dot-dashed line) and the mean temperature Tj^y (dotted 
line) of matter located outside galaxies and clusters. 



ally, at low z the temperature starts decreasing again be- 
cause of the adiabatic expansion of the universe since the 
heating time becomes larger than the Hubble time (see 
upper panel in Fig. 3 below and Sect. 8. 2 in VS). As de- 
scribed in Sect. 2.4 and Valageas et al. (1999a), at low z we 
set the characteristic temperature of Lyman-a clouds to 
3 10'' K when Tjgm becomes smaller than this value. In- 
deed, although the density contrast of small clouds can be 
small and even negative (their "overdensity" (1 -I- A) can 
reach the minimum (1 -I- A)^ ~ 10~^ shown in Fig. 5 be- 
low) they correspond to filaments or underdense objects. 



surrounded by regions of even lower density, which have 
decoupled from the expansion of the universe. Hence they 
do not cool because of adiabatic expansion (or at least this 
term is much smaller than for the IGM). In particular, we 
stress that on strongly non-linear scales even objects with 
a density contrast A equal to are not described by the 
linear theory. Thus, the average density p loses the signif- 
icance it has on large scales in the sense that it does not 
define any longer a density boundary between two physi- 
cally different regimes. In fact, this role is now played by 
the density contrast (1 -I- A)„, defined in (p^), which de- 
scribes most of the volume of the universe seen at small 
scales (see Valageas & Schaeffer 1997). In these regions, 
patches of matter with the density p appear as density 
peaks. This is clearly seen in Valageas et al. (1999b) where 
a comparison with numerical simulations shows that the 
same formulation (^), which is based on a stable-clustering 
approximation in a statistical sense, provides a reasonable 
description of objects defined by a density contrast which 
can vary from (1 -|- A) = 5000 downto (1 -I- A) = 0.5. 

As compared to the original model with only photoion- 
ization heating (upper figure) the main difference when we 
include a source term corresponding to supernovae (mid- 
dle figure) or quasars (lower figure) is that the IGM tem- 
perature keeps increasing until z ~ to reach a value 
Tigm{z = 0) ~ 5 10^ K (the dechne of Tiqm at low z in 
the original case, and for the quasar scenario at z ~ 0, is 
due to adiabatic cooling, because of the expansion of the 
universe). In this case, as seen in (|l2|) we take the Lyman- 
a clouds to be heated to the same temperature in order to 
have a self-consistent model. However, this high tempera- 
ture could make it difficult to recover the observed prop- 
erties of Lyman-a clouds at low z. Although this would 
deserve a detailed study we do not further investigate this 
point in this article where we mainly consider the energy 
requirements implied by efficient reheating. Moreover, the 
opacity due to Lyman-a clouds at low z does not influ- 
ence much the reionization and reheating history of the 
universe (the medium is optically thin as shown by the 
Gunn-Peterson test). Besides, in order to get reliable esti- 
mates of the properties of Lyman-a clouds in the case 
of strong reheating by supernovae or QSOs one would 
certainly need to take into account the spatial inhomo- 
geneities of this reheating process, which is beyond the 
scope of the present study. 

In order to obtain the same entropy and temperature 
at z = for the IGM, we see that the redshift zs (shown 
by the dashed vertical line), where no just-virialized halo 
can cool, is higher for the quasar scenario (QSO) than for 
supernova heating (SN). Indeed, we obtain: 



SN : zs~0.4 



QSO : zs~2 



(36) 



This is due to the peak of the quasar luminosity function 
at z 2 and its sharp decline at lower redshift. Indeed, 
this implies that most of the heating process occurs at 
z ^ 2. On the other hand, since the galaxy luminosity 
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function evolves more slowly and does not drop at 2 ~ 0, 
the heating process due to stars keeps going on until z — 
so that it appears delayed as compared to the quasar sce- 
nario, see discussion in Sect j2.3| . This also shows in the 
behaviour of Tjgm which keeps strongly increasing until 
z = for supernova heating while it remains roughly con- 
stant (and even shows a slight decline) after zs for quasar 
heating. The efhciency factors asN and aq we use are: 



SN : asN = 1 , = 
QSO : asN , ag = 0.008 



(37) 



We can check that they are consistent with the estimates 
(|l|) and (|l|) and the requirement that Tjgm ~ 5 10^ K 
at z = 0. This latter value for Tjqm is due to the con- 
straint {S)iGM — log(A:T?i^^^'^/l keV cm^) ~ 2 which we 
impose in order to get the break at T ^ 1 keV of the rela- 
tion T — Lx for clusters, as explained above and described 
in details in Sect j4.3| . In fact, for supernova heating the ef- 
ficiency factor asAT = 1 is not sufficient (though not by 
far since asN = 1-7 is enough) to explain the break of the 
cluster T—Lx relation, as discussed in Sect.^. However, 
since we must have asAr < 1 we keep this value for the 
supernova efficiency factor. On the other hand, we note 
that a very small value for the quasar efficiency factor ag 
is sufficient to raise the entropy of the IGM to a level high 
enough to explain the cluster observations. Note also that 
aq <C eg, where eg — 0.1 defined in Sect. 2^ measures 
the quasar radiative efficiency. Thus, the quasar scenario 
(QSO) appears to be quite reasonable while the supernova 
hypothesis (SN) seems less likely. However, further work 
is needed in order to assess with a sufficiently good accu- 
racy the efficiency of these two processes of energy transfer 
before one can draw definite conclusions. 



3.2. Entropy production 

We display in Fig.^ the redshift evolution of the entropy 
of the gas. It increases with time as structure formation 
develops and heating processes grow. The specific entropy 
Su of underdense regions is larger than the mean IGM 
entropy {S)igm, and increasingly so at smaller redshift, 
because of their low density rib. As was the case for the 
temperature evolution, we note that the entropy shows a 
faster rise at low z for the supernova scenario (SN) than 
for the quasar heating (QSO). Again this is due to the fact 
that the quasar energy output (proportional to the lumi- 
nosity function) peaks at z ^ 2 contrary to the galaxy con- 
tribution which keeps increasing until z ^ 1. The entropy 
Scooi is defined from the temperature Tcooi and the density 
contrast Agai{z). The fact that it is lower than {S)igm, 
especially at low z, shows that the smallest galaxies are 
influenced by the entropy floor set by the IGM. However, 
we recall that Scooi is not the entropy of the gas in such a 
galaxy, since the gas compression stops when it reaches the 
density contrast Aad < Ac and it subsequently cools and 
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Fig. 2. The redshift evolution of the characteristic en- 
tropies of the universe. We display the mean IGM entropy 
('5') IGM (dot-dashed curve) , the entropy Su of underdense 
regions (solid line) and the entropy which would corre- 
spond to the smallest galaxy Scooi (low dashed line). From 
top downto bottom, the various panels correspond to pho- 
toionization heating only, supernova heating and quasar 
heating, as in Fig.|l|. 

falls into the dark matter potential well, which diminishes 
its entropy (transfered into the radiation). 

We note that if we could measure the temperature and 
the entropy of the gas within Lyman-a clouds or small 
groups at high redshift z ^ 1 we might be able to see which 
scenario, (SN) or (QSO), is best favoured, using the fact 
that the redshift evolution is slower for the quasar heating 
process. However, it is clear that this would not give a 
deflnite answer because of the uncertainty associated to 
these poorly known processes. 

3.3. Time-scales 

We display in FigJ| the redshift evolution of the various 
characteristic time-scales. For the original model (upper 
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Fig. 3. The redshift evolution of the characteristic heat- 
ing and cooling time-scales. All times are shown in units 
of the Hubble time t^f(z). The dashed curve shows the 
cooling or heating time associated to atomic processes 
(Compton cooling, coUisional excitation,...). The solid line 
labelled theat,j shows the photoionization heating time. 
The dot-dashed (resp. dotted) line labelled tE,sN (resp. 
tE,Q) shows the heating time-scale associated to super- 
novae (resp. quasars). From top downto bottom, the var- 
ious panels correspond to photoionization heating only, 
supernova heating and quasar heating, as in Fig.^. 



panel), at large redshift 8 < z < 24 the smallest time- 
scale is the photoionization heating time, defined in (|lj), 
which means that the IGM temperature increases in this 
redshift range (see Figj^). At lower and larger redshift, 
the smallest time is the Hubble time tuiz) (correspond- 
ing to the ordinate in the figure) which implies that the 
IGM cools due to the adiabatic expansion. The dashed 
curve which shows a peak at log(l -\- z) ~ 1.4 corre- 
sponds to atomic processes within the IGM (coUisional 
excitation, coUisional ionization, bremsstrahlung, Comp- 
ton cooling or heating,...). Most of the time it is dominated 
by Compton cooling or heating which explains the peak 



)■ 

At higher redshift the IGM gas is heated by CMB photons 
while at lower z the IGM is cooled through the interac- 
tion with the CMB. Around reionization at z ~ 9 the 
main process is coUisional excitation (see VS for details). 
In both lower panels we also display the heating time-scale 
tE,SN or tE,Q associated to the additional energy source. 
We can see that it becomes the dominant process some- 
what after reionization at z 8. Then, since Ie becomes 
the smallest time-scale the corresponding energy source 
heats the IGM up to T/gm ^ 5 10^ K. We can note again 
that the supernova heating is somewhat delayed as com- 
pared to the quasar scenario. The rise at low z of theat.j 
as compared to the upper panel is due to the growth of 
the IGM temperature, see (|l4|). 

Note that in all three cases the IGM is reheated and 
reionized by the energy output of galaxies and QSOs 
(whether it is radiation or kinetic energy) . Thus in all sce- 
narios we expect a proximity effect (along a line of sight to 
a distant QSO the IGM is more ionized close to the quasar) 
since these processes are not exactly homogeneous (they 
are more efficient close to the sources). This also holds for 
the (SN) scenario because QSOs are associated with galac- 
tic cores (where a black hole is surrounded by an accretion 
disk) which means that they are embedded within regions 
of star formation (which occurs in the host galaxy). In 
fact, this proximity effect provides a measure of the inho- 
mogeneity of the reheating and reionization process. Note 
that Davidsen et al.(1996) and Reimers et al.(1997) do 
not observe any proximity effect which suggests that our 
homogeneous approximation is reasonable. 

3.4- Reionization history 

We show in Fig.^ the redshift evolution of the background 
UV flux J21 defined by: 



J21 



I 10-2^ erg cm-^s-^Hz^^sr-i (38) 



The very sharp rise at z^i — 9 corresponds to the reion- 
ization redshift when the universe suddenly becomes op- 
tically thin. We can see that all scenarios give nearly the 
same results since at large redshift the entropy of the IGM 
is not sufficiently large to significantly affect galaxy and 
quasar formation. This was also apparent in Figj^ and 
Fig.^. In particular, at high z the gas densities are large 
so that cooling is quite easy. As a consequence, the en- 
tropy of the IGM cannot prevent the gas from cooling 
and falling into galactic halos to form stars or quasars so 
that the reionization redshift does not depend on the effi- 
ciency factors a. This is in fact reassuring, since it shows 
that most of the results we obtained in VS (e.g. ioniza- 
tion state of hydrogen and helium,...) are still valid and 
do not depend on the injection of energy into the IGM by 
stars or quasars (the reionization redshift we obtain here 
is larger than in VS because we use Tsn = 10® K instead 
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Fig. 4. Redshift evolution of the background UV flux J21. 
The dotted line corresponds to photoionization heating 
only, the solid line to supernova heating and the dashed 
line to quasar heating. The data points are from Gi- 
allongo et al.(1996) (square), Cooke et al.(1997) (filled 
square), Vogel et al.(1995) (triangle, upper limit), Don- 
ahue et al.(1995) (filled triangle, upper limit) and Kulka- 
rni & FaU(1993) (circle). 



of TsN — 2 10^ K). The UV flux in the quasar scenario 
(QSO) is somevifhat larger than for both other cases at 
2; ~ 7 (log(l + z) ~ 0.9) although the star formation rate 
is a bit smaller (see Fig.^ below) because the opacity of the 
universe is lower. Indeed, this implies a smaller absorption 
of the radiation emitted by stars and quasars. This is due 
to the larger temperature of the IGM (see Figj^) which 
leads to a lower opacity from Lyman-a clouds. However, 
this is only relevant for z ^ 5. 

3.5. Characteristic density contrasts 

We present in Fig.^ the redshift evolution of the character- 
istic density contrasts within the universe. The lower solid 
line shows the "overdensity" (1 -I- A)„ of the underdense 
regions which cover most of the volume, see (p^). At large 
z we have (1 -|- A)tj ~ 1 since very few baryonic structures 
have formed and the universe appears as a nearly uniform 
medium. At low z it declines and can reach very low val- 
ues (H-A)„ ^ 10^^ as "voids" appear amid filaments and 
halos. It is larger for both lower panels because the IGM 
temperature is higher, which means that the scale Rd is 
larger, and dark-matter density fluctuations are smaller 
on larger scales (following the behaviour of the two-point 
correlation function). On the other hand, the mean den- 
sity contrast (1 + A)7G'a/, which takes into account all the 
matter which is not enclosed within galactic halos where 
the gas was able to cool, from "voids" up to filaments and 
forest Lyman-a clouds, remains close to unity. Indeed, the 
fraction of matter 1 — Fcooi which is not embedded within 
galactic halos remains large until z — since we have 
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Fig. 5. The redshift evolution of the characteristic density 
contrasts. We display the density contrast (1 -I- A)^, of 
large underdense regions in the IGM (lower solid line), 
(1 + A)/GM for the mean IGM density (lower dashed line), 
(1 + A)iy for the mass-averaged density within the IGM 
(upper dot-dashed hue) and (1 + A)^^; for galactic halos 
(upper solid line). 



Fcooi Si 0.6 at z = 0, see FigJ|. The mean (1 + A)Ly is 
larger than (1 + A)/G'Af since it is weighted by mass which 
gives more weight to dense regions (filaments, clouds) as 
compared to low-density areas. Finally, the density con- 
trast Agai corresponds to the small galactic halos (which 
are not influenced by the cooling condition R < Rcooi)- 
Thus, at all times in the upper panel and at 2 > 2 in the 
lower panels it is equal to the "virialization" density con- 
trast Ac(z) obtained from the usual spherical model. At 
low redshift in both lower panels it is larger than Ac(z) 
and equal to the expression defined in ( |3^ ) because of en- 
tropy considerations. Again, we note that for the quasar 
scenario (QSO) this effect appears somewhat earlier. 
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Fig. 6. The Compton parameter y up to redshift z de- 
scribing the Sunyaev-Zeldovich effect from the "voids" : yu 
(lower dashed line), the IGM as a whole (voids, filaments 
and clouds): yicM (solid line) and virialized halos: yhaios 
(upper dashed line). 



3.6. Compton y parameter 

The hot gas within the IGM or virialized halos (e.g. clus- 
ters) scatters some photons of the CMB from the low en- 
ergy Rayleigh- Jeans part of the spectrum up to the high 
energy Wien tail. The magnitude of this perturbation is 
conveniently described by the Compton parameter y: 



dt , kT 

C — dz Ut Tie 7 

Q dz rripC'^ 



(39) 



We consider three components to the global Sunyaev- 
Zeldovich effect. Thus we define yu (noted yiGM in VS) 
describing the effect of low-density regions, ywM which 
takes into account all the matter within the IGM (from 
voids up to filaments and Lyman-a forest clouds) and 
Vhaios which describes the effect from virialized halos 
above Tcooi (i-e. galaxies and clusters). Both and yjcM 



reach a plateau at Zri = 9 since at earlier times the uni- 
verse was almost exactly neutral. On the other hand yhaios 
saturates earlier because the fraction of matter embedded 
within massive virialized halos declines at large z. The 
contribution y^ of low-density regions is always small since 
it only contains a small fraction of matter at low z with 
a small temperature. Similarly, yhaios from galaxies and 
clusters is much larger than yicM in the upper panel be- 
cause although both components contain similar fractions 
of baryonic matter the temperature of these virialized ha- 
los is much larger than the IGM temperature (see Fig.|l|). 
However, when there is an additional source of energy from 
supernovae or quasars yicM becomes non- negligible, es- 
pecially for the quasar scenario (QSO) where Tjcm rises 
earlier. Note that on general grounds we can expect yiGM 
to be at most of the same order as yhaios since a large 
fraction of matter is embedded within virialized halos and 
the temperature of the IGM should not be larger than 10^ 
K while the temperature of these collapsed objects can be 
much larger (e.g. clusters have T ~ 10^ K). A more pre- 
cise study of the Compton parameter induced by clusters 
is presented in Valageas & SchaefFer (1999b). In any case, 
we find that the Compton parameter from any contribu- 
tion is much lower than the COBE/FIRAS upper limit 
y < 1.5 lO^'^ provided by observations (Fixsen et al.l996). 
This means that we cannot constrain the entropy scenario 
from its effect on the Compton parameter using current 
observations. On the other hand, it shows that we do not 
contradict the data, hence these heating processes remain 
plausible explanations for the cluster observations. 

We can also compute the X-ray background provided 
by the IGM, from underdense regions up to Lyman-a 
clouds in filaments. However, we find that it is negligi- 
ble since we get at most (for the QSO scenario) a flux 
in the 0.5-2 keV band of $ ~ 10~^'' erg s~^ cm~^ deg~^ 
from the IGM, and $ ^ 10^^^ erg s^^ cm^^ deg^^ from 
virialized halos which have not cooled, as compared to 
the observed extragalactic intensity / ~ 8 10~^^ erg s~^ 
cm~^ deg~^ (Miyaji et al.l998). On the other hand, the 
contribution from galaxies and clusters is of the order of 
10^^^ erg s~^ cm~^ deg~^. Indeed, the X-ray flux is very 
sensitive to the density (it varies as n^) and to the tem- 
perature. Thus, most of the contribution comes from high 
temperature clusters T <; 0.5 keV while the low temper- 
ature of the IGM TiGM 5 10^ K - 0.08 keV impHes a 
very small contribution since cxp(— 0.5/0.08) ^ 10""'^. 

4. Feedback on structure formation 

4-.1. Star formation 

As described in Sect . ^.5. 3 the entropy of the IGM inhibits 
the cooling of the gas which in turns decreases the effi- 
ciency of star formation. In particular, after the redshift 
zs the gas which would become embedded within new non- 
linear structures cannot cool within a few Hubble times at 
formation because its entropy prevents its density to reach 
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large enough values to start cooling efficiently. Thus, in 
our model we obtain a population of "old" galaxies which 
gradually convert their matter content into stars but the 
overall fraction of gas which can cool does not increase any 
more. After a while, when a large part of this matter has 
formed stars, this leads to a decrease of the star formation 
rate at low z. Note that we neglect cooling flows within 
groups and clusters which provide an additional source of 
cool gas which may form stars. However, this feedback ef- 
fect onto galaxy formation is likely to persist in a more 
detailed model. 
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— 1 — 



-2 
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Fig. 7. The redshift evolution of the comoving star forma- 
tion rate. The dotted line corresponds to photoionization 
heating only, the solid line to supernova heating (SN) and 
the dashed line to quasar heating (QSO). The data points 
are from Madau (1999), see references therein. 

We show in Fig j?] the comoving star formation rate we 
obtain for the three scenarios. We clearly see the inhibi- 
tion of star formation as compared to the case with pho- 
toionization heating only. This also shows at large redshift 
z > zg, where new galaxies appear but the high entropy 
of the gas prevents the formation of the smallest galaxies 
which occur in the case with (Xsn = olq = 0. As ex- 



plained in Sect. 2.5.3 this feedback effect is larger at low 
z where structure formation is further developped (hence 
the energy source is high) and the baryonic density is lower 
(hence cooling is less efficient). Moreover, this feedback ef- 
fect onto star formation may partly explain the sharp drop 
of the star formation rate observed at low redshift z <\. 
Indeed, this decline is not very easy to obtain in usual 
models since at z ~ the fraction of baryons which has 
been converted into stars is still very small F^tar ^ 10%. 
As a consequence, it is difficult to get a sudden stop of 
the star formation process because there is plenty of gas 
available (although one can obtain in a natural fashion a 
decline as shown by the dotted line). On the other hand, 
the supernova or quasar heating of the IGM is able to 
suddenly change the conditions of star formation when 
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Fig. 8. The redshift evolution of the fraction of baryonic 
matter embedded within cooled objects Fcooi (solid line) 
and stars Fstar (dashed line). The data point at z ~ 
shows the observed mass within stars, from Fukugita et 
aL(1998). 



the entropy of the gas becomes of the order of the en- 
tropy generated by gravitational collapse within the new 
non- linear halos. 

Although the shape of the redshift evolution of the 
comoving star formation rate we obtain for the quasar 
scenario (dashed line) agrees with observations, its nor- 
malization is somewhat too low. This might be the sign 
of a shortcoming of our description, in particular the as- 
sumption that reheating is uniform could lead to an over- 
estimate of the entropy feedback onto star formation, since 
in a more realistic model with inhomogeneous reheating 
this effect may require more time to affect all the gas (es- 
pecially for the small galaxies which form far away from 
clusters or groups of large galaxies and quasars). In fact, 
the star formation rate for the photoionization only case 
is already a bit two low at z ^ 1. Thus our model of 
star formation is not perfect yet and we might underesti- 
mate the supernovae heating by a factor 2. Although this 
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would translate into a supernovae efRciency factor smaller 
than unity, asN — 0.8, this value remains quite large and 
it does not modify our conclusions (e.g. that the IGM is 
more likely to have been reheated by quasars). Thus, we 
think these results are already quite encouraging, in view 
of the simplicity of our model. On the other hand, note 
that we obtained a correct amplitude for the UV flux as 
shown in Fig.^. 

We show in Fig J| the fraction of matter enclosed within 
cooled objects (galaxies) and stars {Fstar — ^star/^b)- We 
see that our results agree with the observed mass of stars 
(which may be more robust than the observed redshift evo- 
lution of the star formation rate itself). This could be ex- 
pected since our model for galaxy formation (though in a 
more detailed version) was already checked in Valageas & 
Schaeffer (1999a) against observations. In particular, the 
cooled fraction Fcooi sd z — obtained with photoioniza- 
tion heating only (upper panel) is rather large. However, 
the mass of stars we obtain agrees with observations which 
shows that we do not encounter the overcooling problem. 
This is due to the small star formation efficiency implied 
by several processes: ejection of matter by supernovae and 
stellar winds, energy released by halos mergings and col- 
lapse. 

Of course, we recover the behaviour seen in Fig.^. In 
particular, for strong heating of the IGM we clearly see 
the saturation of the mass of cooled gas at low z and the 
upper bound for the mass of stars it implies. Note that 
this feedback effect is quite general. Indeed, we require 
the reheating of the IGM to be large enough to affect 
cluster formation at z ~ for halos with a virial temper- 
ature T ^ 0.5 keV. It is clear that this implies a signif- 
icant effect onto galaxy formation at low redshift, since 
galaxies consist of shallower potential wells T ^ 0.1 keV. 
Since we must simultaneously describe galaxies, quasars 
and clusters, this leads to contradictory constraints. In- 
deed, the cluster T — Lx relation requires a high reheat- 
ing temperature T ^ 5 10^ K (see Sect. 4.3) while the 
observed star formation rate requires a small enough re- 
heating, T ^ 5 10^ K, so as not to inhibit too much galaxy 
formation. Thus, we see from Fig.^ and Fig.^that these 
constraints imply a reheating temperature T ~ 5 10^ K 
and ag ~ 0.008. This clearly shows the importance of us- 
ing global models like ours (even though simplified) which 
allow one to evaluate the consequences of such processes 
on all objects (galaxies, clusters,...). Indeed, this provides 
strong constraints on such descriptions and it is the only 
way to test the global validity of these scenarios which 
should simultaneously account for all structure formation 
processes. We also note that a large fraction of baryonic 
matter at z ~ is embedded within density fluctuations 
in the IGM and Lyman-a clouds with a rather large tem- 
perature T ^ 5 10^ K. Part of this component may be 
difficult to observe. 



4-2. Quasar luminosity function 




Fig. 9. The evolution with redshift of the B-band quasar 
luminosity function in comoving Mpc^3. The dotted hues 
correspond to photoionization heating only, the solid lines 
to supernova heating and the dashed lines to quasar heat- 
ing. The data points are from Pei (1995). 

We display in Fig.^ the redshift evolution of the quasar 
luminosity function for all scenarios. We can see that at 
large redshifts z ^ 2 all curves nearly superpose since 
the additional energy source (from quasars or supernovae) 
plays no role at early times. On the other hand, as was the 
case for star formation, at small redshift z <, 0.5 the en- 
tropy "floor" induced by quasar heating leads to a decrease 
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of the quasar luminosity function which slightly improves 
the agreement with observations. Thus, although the de- 
cline at low z of the QSO multiplicity function is a natural 
outcome of our model even for the original scenario with 
photoionization heating only (dotted line), part of this 
decrease may also be due to the entropy production by 
quasars at earlier times. In view of the simplicity of our 
model for quasar formation, which is a natural outcome of 
our description of galaxies, we think our predictions agree 
reasonably well with observations. Of course, the physics 
of quasars may be more intricate than the description we 
use in this article but any meaningfuU improvement would 
require a detailed model of the accretion processes lead- 
ing to quasar formation, which is beyond the scope of this 
study. 

J^.S. Clusters 

As we noticed in the introduction, the reheating of the 
IGM by supernovae or quasars at z J; 1 can affect the 
formation of clusters at lower redshifts z ^ 0.5 since it 
leads to a mimimum entropy {S)igm of the gas which can 
break the expected scaling of the relation temperature - X- 
ray luminosity of clusters. In order to obtain an estimate 
of the characteristic virial temperature Tad.ci where this 
transition should occur we define: 

Tad.ci = TJGM (}^) 10<^)-"-^" (40) 

where we used (^2|) . This is the temperature of the gas at 
the density contrast Ac(z) (which defines clusters in our 
model) with the mean entropy of the IGM {S)igm- Thus, 
massive clusters with T ^ Tad.ci are not affected by the 
entropy floor of the IGM because shock heating during the 
gravitational collapse of the halo generates a much larger 
entropy so that the usual scaling law is recovered. On the 
other hand, within smaller potential wells the gas has a 
larger entropy than the one produced by shock heating 
which leads to a smoother gas density profile and to a 
smaller density. This implies a smaller luminosity since 
Lx oc rig . In (|o|) we used the virial density rather than the 
core density, which is significantly larger, because most of 
the mass is characterized by densities of order Ac and we 
assume isothermal equilibrium. However, it is clear that 
(|^) is only a rough estimate which could be uncertain by 
a factor 2. 

We show in Fig.|l^ the redshift evolution of Tad,ci- 
Again, we can check that in the quasar scenario (QSO) the 
heating of the IGM occurs earlier than for the supernova 
case (SN). Moreover, since the quasar luminosity function 
drops at low z so that adiabatic cooling becomes the main 
process, as seen in Fig.^, Tad,ci decreases at low z. On the 
other hand, for the supernova scenario the fact that the 
redshift zs is smaller leads to a smoother evolution at 
z '-^ 0. In any case, we see that we obtain at small red- 
shift a characteristic temperature Tad.ci ~ 0.5 keV. This 
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Fig. 10. Redshift evolution of the characteristic temper- 
ature Tad,ci which describes the influence of reheating on 
cluster formation. The dotted line corresponds to pho- 
toionization heating only, the solid line to supernova heat- 
ing and the dashed line to quasar heating. 

is similar to the values which are used in studies of cluster 
formation (Cavaliere et al.l997 use Tad,ci — 0.5 — 0.8 keV 
while Valageas & Schaeffer 1999b use Tad.ci — 0.35 keV). 
Indeed, this is how we chose the parameters a. Of course, 
smaller efficiency factors than those set in (^^ lead to 
smaller IGM entropy and lower Tad,ci- However, it is in- 
teresting to note that a larger aq does not change much 
our results at z = for the quasar scenario (QSO) be- 
cause of the feedback effect we described above. On the 
other hand, for the supernova scenario (SN) the available 
range is already limited by the upper bound asN < 1- 

From the mean entropy of the IGM, or the characteris- 
tic temperature Tad,ci, we can obtain the cluster tempera- 
ture - X-ray luminosity relation, as in Valageas & Schaeffer 
(1999b). The bolometric X-ray luminosity Lx of a cluster 
of volume V is: 

Uoi = I nlA,iTg)dV (41) 
Jv 

where rig is the electron number density and Ah{Tg) is the 
bremsstrahlung emissivity function (in erg cm'^ s~^) at 
temperature Tg. Thus, contrary to the Sunyaev-Zel'dovich 
effect, see (|9|), the X-ray luminosity strongly depends on 
the density profile of the hot gas within the cluster. During 
the gravitational collapse of the cluster shocks heat the 
gas up to the virial temperature T of the dark matter 
halo. However, the adiabatic compression of the gas from 
the IGM also heats the gas up to Tad,ci which is a lower 
bound for the gas temperature Tg. In order to take into 
account both of these effects, we simply write for the final 
temperature of the gas: 

Tg = T + Tad,ci (42) 

Next, in order to obtain the density profile of the gas 
within the dark matter potential well we assume isother- 
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mal and hydrostatic equilibrium at the temperature Tg 
and we obtain (see also Cavaliere & Fusco-Femiano 1978): 



Pg oc p oc r 



with 



T 



T 



Tg T + Tad,c 



(43) 



where we used an isothermal profile p oc r^^ for the dark 
matter halo. Thus, for deep potential wells T ^ Tad,ci we 
have /3 ~ 1 and the gas follows the dark matter density 
profile while for cool clusters T ^ Tad we get (3^1 and 
the gas density profile is smoother than the dark matter 
distribution. This change of the shape of the gas density 
profile leads to a break in the relation Tg — Lx- More- 
over, in the inner parts of the cluster the density is large 
enough to lead to a small cooling time so that a cooling 
flow develops and some of the gas forms a cold component 
which does not emit in X-ray any longer. Thus, we define 
the cooling radius Rc as the point where the gas density 
reaches the threshold pgc such that: 



tcooi = tn with tcooi — 



iplnipkTg 

2pPgcKc{Tg) 



(44) 



where Kc{Tg) is the cooling function (which is dominated 
by bremsstrahlung cooling for T > 1 keV) and tniz) is 
the Hubble time. At large radii r > Rc the density is 
lower than pgc hence the local cooling time is larger than 
the Hubble time. Then, the gas distribution and the tem- 
perature had not had time to evolve much and the X-ray 
emissivity is proportional to p"^ Ab{Tg), see (p|). On the 
other hand, within the cooling radius Rc the gas had time 
to cool and form dense cold clouds. However, we consider 
that some of the gas is still hot and emits in X-ray as cool- 
ing does not proceed in a uniform fashion (Nulsen 1986; 
Teyssier et al.l997). The density of this remaining gas has 
to be of order pgc and we write the X-ray luminosity of 
the cluster as: 



Lx = e 



1 - 




4/3-3' 



(45) 



where the factor 1 describes the contribution of the core, 
within Rc, and the second term comes from the halo (note 
that the contribution within Rc is never much larger than 
the one from r ^ Rc). The factor e = 3 is a parameter 
of order unity which we use to normalize our relation to 
observations for massive clusters {Tg > 1 keV). We ex- 
pect e ^ 1 which is indeed the case. It measures the den- 
sity fluctuations of the gas distribution, since Lx oc (rig) 



and at any radius 



> 



\ Note that our descrip- 



tion is similar to the model developed by Cavaliere et 
al. (1997, 1998) to describe the relation between the gas 
and the dark matter density profile, which shows in the 
quantity /3. However, while they define a core radius from 
the density distribution itself (because they use a dark 



matter profile which grows more slowly than at small 
r) the core radius we use describes the cooling of the gas, 
independently of the profile of the underlying dark matter 
halo. Using a shallower dark matter density profile would 
give similar results with a slightly larger Tad,ci and aq 
(see also Valageas & Schaeffer 1999b). 
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Fig. 11. The cluster temperature - X-ray luminosity rela- 
tion at z = (solid lines) and z — 1 (dashed lines). The 
upper panel shows the supernova scenario (SN) and the 
lower panel the quasar scenario (QSO). The data points 
are from Mushotzky & Scharf (1997) for clusters and from 
Ponman et al.(1996) for groups. 

We show our results in Figjl^ at z = and z = I, for 
both (SN) and (QSO) scenarios, using the redshift evolu- 
tion we obtained in Fig.|l^ for Tad,ci- First, we can check 
that our results agree with observations for hot clusters 
T > 1 keV. Then, we see that the initial entropy of the 
IGM leads to a break of the relation Tg — Lx at low tem- 
peratures. However, for the supernova heating scenario we 
would require asN — 1-7 to get a sufficiently large effect 
(though even for agA? = 1 we already see a knee in the re- 
lation T — Lx)- On the other hand, for the quasar heating 
case we need aq ~ 0.008. Note that the strong redshift 
evolution of Tad.ci in the quasar scenario, seen in Fig.p^ 
leads to a clear redshift dependence of the break of the 
relation Tg — Lx ■ On the contrary, the slow redshift evolu- 
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tion of Tad,ci in the supernova scenario leads to a smoother 
redshift dependence of the relation Tg — Lx ■ This suggests 
that observations of the evolution of the temperature - lu- 
minosity relation at low Lx could provide some contraints 
on the (QSO) scenario. On the other hand, at large Lx 
there is almost no redshift evolution, which is consistent 
with observations (Mushotzky & Scharf 1997). Finally, we 
note that we may underestimate the effect of the heat- 
ing from supernovae or quasars onto cluster formation in 
our model. Indeed, while we assumed this energy source 
to be homogeneous the gas which will later build a clus- 
ter is more likely to be reheated than an average calcu- 
lation would show since clusters form at density peaks 
where the local density of galaxies and quasars is higher 
than average. This means that the break of the relation 
Tg — Lx could appear at a slightly larger temperature than 
shown in Fig.|l^. This might "help" the supernova sce- 
nario as an efficiency factor asN smaller than unity could 
be sufficient. However, it would probably remain close to 
(XsN = 1- On the other hand, we note that our constraints 
for the reheating of the IGM do not depend much on the 
details of the model of clusters since in any case in order 
to get a break of the Tg — Lx relation at Tg ^ 1 keV one 
necessarily needs to introduce a characteristic tempera- 
ture Tad,ci ~ 0.5 keV (FigJlO|) which sets the location of 
this bend. A more detailed model of clusters is presented 
in Valageas & Schaeffer (1999b) where a good match with 
observations is obtained with Tad,c — 0.35 keV. 

5. Conclusion 

In this article we have described an analytic model for 
structure formation processes which deals in a consistent 
fashion with quasars, galaxies, Lyman-a absorbers and 
underdense regions within the IGM. This allows us to ob- 
tain the reheating and reionization history of the universe, 
as well as the evolution of the entropy of the gas. We 
considered three scenarios, with different efficiency fac- 
tors for the transfer of energy from supernovae or quasars 
into the IGM. Thus, we have shown that the energy pro- 
vided by quasars is sufficient to reheat the universe and 
raise the mean entropy of the IGM up to the value re- 
quired to match the "floor" level observed in cool clus- 
ters. This is an upper bound on the entropy production 
and this value allows to explain the behaviour of the clus- 
ter T — Lx relation. On the other hand, the supernova 
heating scenario would require an efficiency factor of or- 
der unity {asN — 1-7). Thus, the IGM is more likely to 
have been reheated by quasars than by supernovae. 

Of course, a more realistic treatment would account for 
the details of the quasar interaction with their gaseous en- 
vironment. However, our study already shows that quasar- 
driven outflows can provide an important heating mecha- 
nism. On the other hand, a detailed model of the inhomo- 
geneous character of this reheating is probably necessary 
in order to evaluate its effects on Lyman-a clouds since 



most of the opacity may come from clouds located far 
away from quasars which have not been reheated. 

We showed that the feedback of entropy production 
onto structure formation may partly account for the de- 
cline at low z of the comoving star formation rate and of 
the quasar luminosity function, in addition to cluster ob- 
servations. This is an interesting prospect since it links 
different processes to the same phenomenon and it gives 
additional weight to the hypothesis of such a reheating 
scenario. Moreover, it provides a narrow range for the re- 
heating of the IGM (T ~ 5 10^ K) since we must satisfy 
the contradictory constraints arising from clusters (which 
require a large reheating so as to modify the T — Lx re- 
lation) and from galaxies and quasars (which require a 
small reheating so that galaxy formation is not too much 
inhibited). Thus, it is important to simultaneously address 
these processes in order to check the validity of a given sce- 
nario. On the other hand, the reionization process of the 
universe is almost not modified which means that our re- 
sults for the ionization state of the gas and the background 
UV flux should be quite robust. 

Then, we found that although both scenarios, normal- 
ized to the current value of the entropy measured in cool 
clusters, are very similar, the reheating due to quasars oc- 
curs a bit earlier {zs ~ 2) than for supernovae (zs — 0.4) 
because of the sharp drop at z < 2 of the quasar luminos- 
ity function. This might give a clue to discriminate these 
two sources of energy from observations. However, it is 
clear that further work is needed in order to get stronger 
constraints on the possible efficiency of these reheating 
processes, for instance through very detailed numerical 
simulations. Nevertheless, despite the small discrepancies 
we get for the star formation rate or the quasar luminosity 
function as compared with observations, it appears quite 
remarkable that a simple analytic model like ours, which 
necessarily involves some approximations (e.g. we do not 
take into account the scatter of the galaxy or quasar prop- 
erties, nor the inhomogeneity of the supernova or quasar 
heating), provides such a good description of structure 
formation processes. Indeed, we note that at z ~ we 
describe objects which span a wide range in density, from 
(1-|- A) ~ 10"^ for voids and low-column density Lyman-a 
absorbers up to (l-f A) ~ 10'^ for old galaxies, as well as in 
mass, from lO^Af0 for Lyman-a clouds up to IO^^Mq for 
clusters. The fact that we can build a unified consistent 
model for this broad variety of structures strongly suggests 
that hierarchical scenarios like ours, with adequate mod- 
els for galaxy formation and radiative processes, provide 
a realistic description of the actual universe. 
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